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\ We reconsider neutrino decay as an explanation for atmospheric neutrino 

Q^ ■ observations. We show that if the mass-difference relevant to the two mixed 

states and v T is very small (< 10~ 4 eV 2 ), then a very good fit to the 
observations can be obtained with decay of a component of Vn to a sterile 
Qh! neutrino and a Majoron. We discuss how the K2K and MINOS long-baseline 

\ experiments can distinguish the decay and oscillation scenarios. 
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Recently a neutrino decay explanation for the atmospheric neutrino observations of 
Super-Kamiokande was proposed [0. To recapitulate briefly, in presence of a decaying 
neutrino state, the survival probability of is given by 

Pfa -> u„) = sin 4 9 + cos 4 6e~ aL / E + 2 sin 2 9 cos 2 9e~ aL ^ E cos (^§^) (1) 

where we have considered only two state mixing for simplicity: 8m\^ = m\ — m 2 , a = 1712/T2 
and Up = cos9v2 + sin Qv$. If the decay of V2 is z^2 — ^ ^3 + J, where J is a massless scalar, 
then the Sm 2 in the decay is the same as in Eq. (1) and it can be shown that this Sm 2 has 
to be larger than 0.73 eV 2 to satisfy constraints from K — > fi + neutrals decay 0. Then the 
oscillating term in Eq. (1) averages to zero and Piv^, — > v^) simplifies to 

P(v» -> Vy) = sin 4 9 + cos 4 9e~ aL/E (2) 

This decay scenario was analyzed in Ref. [0J and fit to the L/E u dependence and the asym- 
metry of the contained events in Super-Kamiokande [[|. A satisfactory fit with cos 2 9 ~ 0.87 
and a ~ 1 GeV / (where -De is the diameter of the earth) was found. However, it has since 
been shown that the fit to the higher energy events in Super-K (especially the upcoming 
muons) is quite poor ||, 0. 

Another possibility, mentioned in Ref. |L), is that the decay of v 2 is into a state with 
which it does not mix. For example, the three weak coupling states u^, u T , v s (where v s is a 
sterile neutrino) may be related to the mass eigenstates ^2, V3, v\ by the approximate mixing 

/ v {t \ ( cos ^ sin ^ \ / z/ 2 \ 
v T \ = \ -sin^ costf 1/3 (3) 

V Ps J V 1 / V ^4 / 

and the decay is v-i — > z/4 + J. The electron neutrino, which we identify with u\, cannot mix 
very much with the other three because of the more stringent bounds on its couplings [0, 
and thus our preferred solution for solar neutrinos would be small angle matter oscillations. 

Then the 5m| 3 in Eq. (1) is not related to the 5m\ A in the decay, and can be very small, 
say < 10 _4 eV 2 (to ensure that oscillations play no role in the atmospheric neutrinos). In 
that case, the oscillating term is 1 and Piv^ —> v^) becomes 

_> Uft ) = (sin 2 9 + cos 2 9e~ aL ' 2E ) 2 (4) 

This is identical to Eq. (13) in Ref. JTJ. Here we consider this decay model and compare it 
to observations. 

In order to compare the predictions of this model with the standard <-> v T oscillation 
model, we have calculated with Monte Carlo methods the event rates for contained, semi- 
contained and upward-going (passing and stopping) muons in the Super-K detector, in the 
absence of 'new physics', and modifying the muon neutrino flux according to the decay or 
oscillation probabilities discussed above. We have then compared our predictions with the 
SuperK data ||, calculating a \ 2 t° quantify the agreement (or disagreement) between data 
and calculations. In performing our fit (see Ref. || for details) we do not take into account 
any systematic uncertainty, but we allow the absolute flux normalization to vary as a free 
parameter (3. 
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The 'no new physics model' gives a very poor fit to the data with \ 2 — 281 for 34 d.o.f. 
(35 bins and one free parameter, (3). For the standard <-> v T oscillation scenario the best fit 
has x 2 — 33.3 (32 d.o.f.) and the values of the relevant parameters are Am 2 = 3.2 x 10~ 3 eV 2 , 
sin 2 29 = 1 and (3 = 1.15. This result is in good agreement with the detailed fit performed by 
the SuperK collaboration || giving us confidence that our simplified treatment of detector 
acceptances and systematic uncertainties is reasonable. The decay model of Equations (3) 
and (4) above gives an equally good fit with a minimum \ 2 = 33.7 (32 d.o.f.) for the choice 
of parameters 

r u /m u = 63 km/GeV, cos 2 9 = 0.30 (5) 

and normalization (3 = 1.17. 

In Fig. 1 we compare the best fits of the two models considered (oscillations and decay) 
with the SuperK data. In the figure we show (as data points with statistical error bars) the 
ratios between the SuperK data and the Monte Carlo predictions calculated in the absence 
of oscillations or other form of 'new physics' beyond the standard model. In the six panels we 
show separately the data on e-like and /x-like events in the sub-GeV and multi-GeV samples, 
and on stopping and passing upward-going muon events. The solid (dashed) histograms 
correspond to the best fits for the decay model (y^ <-> v T oscillations). One can see that the 
best fits of the two models are of comparable quality. The reason for the similarity of the 
results obtained in the two models can be understood by looking at Fig. 2, where we show 
the survival probability Piy^ — > z/ M ) of muon neutrinos as a function of L/E v for the two 
models using the best fit parameters. In the case of the neutrino decay model (thick curve) 
the probability Piy^ — > v^) monotonically decreases from unity to an asymptotic value 
sin 4 # ~ 0.49. In the case of oscillations the probability has a sinusoidal behaviour in LjE v . 
The two functional forms seem very different; however, taking into account the resolution 
in L/E u , the two forms are hardly distinguishable. In fact, in the large L/E v region, the 
oscillations are averaged out and the survival probability there can be well approximated 
with 0.5 (for maximal mixing). In the region of small L/E„ both probabilities approach 
unity. In the region L/E u around 400 km/GeV, where the probability for the neutrino 
oscillation model has the first minimum, the two curves are most easily distinguishable, at 
least in principle. 

Decay Model 

There are two decay possibilities that can be considered: (a) v 2 decays to P4 which is 
dominantly v s with v 2 and u 3 mixtures of and u T , as in Eq. ([3|), and (b) v 2 decays into 
z/4 which is dominantly v T and v 2 and u 3 are mixtures of and v s . In both cases the decay 
interaction has to be of the form 

£int = #24^" V 1 L J + k.C. (6) 

where J is a Majoron field that is dominantly iso-singlet (this avoids any conflict with the 
invisible width of the Z) . Viable models for both the above cases can be constructed ||,[7| . 
However, case (b) needs additional iso-triplet light scalars which cause potential problems 
with Big Bang Nucleosynthesis (BBN), and there is some preliminary evidence from SuperK 
against z/^-za, mixing 0. Hence we only consider case (a), i.e. u 2 — > v± + J with z/ 4 w u s , as 
implicit in Eq. (|3|). With this interaction, the v 2 rest-lifetime is given by 
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= 16tt m 2 , . 

T2 # 2 ' 5m 2 (l + x) 2 ' 1 ' 

where 5m 2 = m 2 — m\ and x = 7714/1712 (0 < x < 1). From the value of a^ 1 = T2/7712 = 
63 km/GeV found in the fit and for x = 0, we have 

g 2 5m 2 ~ 0.16 eV 2 (8) 

Combining this with the bound on g 2 from K /j, decays of g 2 < 2.4 x 10~ 4 || we have 

5m 2 > 650 eV 2 . (9) 

Even with a generous interpretation of the uncertainties in the fit, this Sm 2 implies a min- 
imum mass difference in the range of about 25 eV. Then v 2 and 1/3 are nearly degenerate 
with masses ^ (9(25 eV) and is relatively light. We assume that a similar coupling of 
v-& to z/4 and J is somewhat weaker leading to a significantly longer lifetime for z/3, and the 
instability of z/3 is irrelevant for the analysis of the atmospheric neutrino data. 

For the atmospheric neutrinos in SuperK, two kinds of tests have been proposed to 
distinguish between v il -v T oscillations and v^-Vg oscillations. One is based on the fact that 
matter effects are present for v^—Vg oscillations |J but are nearly absent for Vy,—v T oscillations 
Toll leading to differences in the zenith angle distributions due to matter effects on upgoing 



neutrinos ||11|| . The other is the fact that the neutral current rate will be affected in v^-v s 



oscillations but not for — v T oscillations as can be measured in events with single 7r s |T2 



In these tests our decay scenario will behave as a hybrid in that there is no matter effect 
but there is some effect in neutral current rates. 

Long-Baseline Experiments 

The survival probability of z/^ as a function of L/E is given in Eq. (1). The conversion 
probability into v T is given by 

P( U(t v T ) = sin 2 6 cos 2 9(1 - e ~ aL/2E ) 2 . (10) 

This result differs from 1 — Piv^ —>■ v^) and hence is different from v^-Pr oscillations. 
Furthermore, Piv^ —+ v^) + P(v^ — > u T ) is not 1 but is given by 

- + P(^ -+ Vr ) = l- cos 2 9(1 - e- aL ' E ) (11) 

and determines the amount by which the predicted neutral-current rates are affected com- 
pared to the no oscillations (or the v^-Vr oscillations) case. In Fig. 3 we give the results for 
P{ u n ~> v tj)i P{ v fi u t) and P(vn — > v^) + P(Vfi — > u T ) for the decay model and compare 
them to the v^-v? oscillations, for both the K2K |13| and MINOS 0] (or the corresponding 



European project |15|) long-baseline experiments, with the oscillation and decay parameters 
as determined in the fits above. 

The K2K experiment, already underway, has a low energy beam E v w 1-2 GeV and a 
baseline L = 250 km. The MINOS experiment will have 3 different beams, with average 
energies E v = 3, 6 and 12 GeV and a baseline L = 732 km. The approximate L/E u ranges 
are thus 125-250 km/GeV for K2K and 50-250 km/GeV for MINOS. The comparisons in 
Figure 3 show that the energy dependence of survival probability and the neutral current 
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rate can both distinguish between the decay and the oscillation models. MINOS and the 
European project may also have r detection capabilities that would allow additional tests. 

Big Bang Nucleosynthesis 

The decay of v-i is sufficiently fast that all the neutrinos (u e , z/^, z/ T , v s ) and the Majoron 
may be expected to equilibrate in the early universe before the primordial neutrinos decouple. 
When they achieve thermal equilibrium each Major ana neutrino contributes N v = 1 and the 
Majoron contributes N v = 4/7 giving and effective number of light neutrinos N v = 4| 
at the time of Big Bang Nucleosynthesis. From the observed primordial abundances of 4 He 
and 6 Li, upper limits on N v are inferred, but these depend on which data are used |T7|-|r9|. 
Conservatively, the upper limit to N v could extend up to 5.3 (or even to 6 if 7 Li is depleted 
in halo stars AT/ 



Cosmic Neutrino Fluxes 

Since we expect both v-i and z/3 to decay, neutrino beams from distant sources (such as 
Supernovae, active galactic nuclei and gamma-ray bursters) should contain only v e and v e 
but no Ufj,, Ufj,, v T and u T . This is a very strong prediction of our decay scenario. 

Reactor and Accelerator Limits 

The u e is essentially decoupled from the decay state z/ 2 so the null observations from the 
CHOOZ reactor are satisfied ||20|| . The mixings of and v T with u s and u e are very small, so 



there is no conflict with stringent accelerator limits on flavor oscillations with large 5m [21 
Conclusions 

In summary, we have shown that neutrino decay remains a viable alternative to neutrino 
oscillations as an explanation of the atmospheric neutrino anomaly. The model consists of 
two nearly degenerate mass eigenstates z/2, z/3 with mass separation > (9(25 eV) from another 
nearly degenerate pair 1/1, z/4. The and z/ r flavors are approximately composed of z/2 and 
z/3, with a mixing angle 6*23 — 57°. The state z/ 2 is unstable, decaying to z/4 and a Majoron 
with a lifetime r 2 ~ 10~ 12 sec. The electron neutrino z/ e and a sterile neutrino u s have 
negligible mixing with z/ M , v T and are approximate mass eigenstates (z/ e W!/i,z/ s « z/4), with 
a small mixing angle #14 and a 5m\ x ~ 10 -5 eV 2 to explain the solar neutrino anomaly. The 
states z/3 and z/4 are also unstable, but with z/3 lifetime somewhat longer and z/4 lifetime much 
longer than the z/ 2 lifetime. This decay scenario is difficult to distinguish from oscillations 
because of the smearing in both L and E v in atmospheric neutrino events. However, long- 
baseline experiments, where L is fixed, should be able to establish whether the dependence 
of L/E v is exponential or sinusoidal. In our scenario only v\ is stable. Thus, neutrinos of 
supernovae or of extra galactic origin would be almost entirely v e . The contribution of the 
electron neutrinos and the Majorons to the cosmological energy density Q is negligible and 
not relevant for large scale structure formation. 
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FIG. 1. Comparison of decay model (solid histograms) and v^-Vt oscillation model (dashed 
histograms) with SuperK data from Ref. ||. 



7 




12 3 4 



x = Log 10 [L/E (km/GeV)] 



FIG. 2. Survival probabiliity for versus log 10 (L / E) for the decay model (heavy solid curve) 
and oscillation model (thin curve). 
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FIG. 3. Long-baseline expectations for the K2K and MINOS long-baseline experiments from 
the decay model and the v^-v-t oscillation model. The upper panel gives the neutral current 
predictions compared to no oscillations (or v^—Vr oscillations). 
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